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We study a three-dimensional abelian lattice model in which the analogue of a theta term 
J>~Y can be defined. This term is defined by introducing a neutral scalar field and its effect 
is to couple magnetic monopoles to the scalar field and vortices to the gauge field. An 
^ | interesting feature of this model is the presence of an exact duality symmetry that acts on 
a three parameter space. It is shown that this model has an interesting phase structure 
for non-zero values of 6. In addition to the usual confinement and vortex phases there are 

O ' 

phases in which loops with composite charges condense. The presence of novel point like 
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excitations also alters the physical properties of the system. 
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1 Introduction 



Topology has played an important role in statistical mechanics and field theory. The 
topology of the dynamical variables leads to novel excitations that can have a profound 
effect in determining the physical properties of the system. A well known example in 
statistical mechanics is the vortex excitation [|TJ in the two-dimensional planar model. These 
vortex excitations, which exist because of the angular nature of the spin variables, drive a 
phase transition which is very different from other phase transitions in statistical mechanics. 
In gauge theories the analog of the vortex excitation is the magnetic monopole ||. The 
magnetic monopole exists as a consequence of the topology of the gauge group. Mult- 
monopole states are also present as collective excitations. These multi-monopole excitations 
can form a plasma phase in which electric charges are confined. A classic example of 
this phenomenon is compact QED in three Euclidean dimensions ||. Apart from the 
intrinsic topology of the dynamical variables, we can also consider terms which have a 
direct topological significance. One such term is the well known 9 term in non-abelian 
gauge theories. For the 577(2) non-abelian gauge theory this term is defined as 

J d 4 x tr F^F^ . (1) 

Although this term is a total derivative it can have a non-trivial effect whenever long range 
fields are present. In the space of finite action configurations, the above term is precisely 
the winding number of mappings from S3 to S3. Gauge field configurations which give a 
non-zero value for the expression in Eq. [1] can affect some of the physical properties of the 
system M. An important effect of the term in Eq. [I] is that it can convert a magnetic 
monopole into a particle with an additional electric charge - a dyon ||. It was conjectured 
in || that non-abelian gauge theories could have new phases of oblique confinement as 
a result of the interactions between these dyons. The analog of the term in Eq. [l| in an 
abelian gauge theory is 

J d 4 xF^ v . (2) 
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Unlike as in the non-abelian theory, this term does not have a topological significance as 
the winding number of any mapping. However, this term can have a non-trivial physical 
effect in the presence of magnetic monopoles. The effect of the 9 term in Eq. |2| was studied 
on the lattice in 0. It was shown that the 9 term drastically alters the phase structure 
of the theory and a rich phase structure was uncovered as a function of 9. The oblique 
confinement phases conjectured in || were also elucidated. An exact duality symmetry ( 
which is the action of the group SL(2, Z) ) was demonstrated to hold in this model || and 
the action of this symmetry was used to predict the entire phase structure of the model. 
Both the 9 terms discussed so far require four Euclidean dimensions. In |7|, § some two 
dimensional spin models were also considered in which a 9 term could be defined. We 
would like to consider such a term in three Euclidean dimensions and study its effects. An 
inspection of the properties of the e tensor shows that defining gauge invariant terms using 
only pure gauge fields is not possible. An exception is the Cherns-Simons term about which 
we will comment later. The 9 term that we will consider in three dimensions is the lattice 
analogue of the following term 

i9 J d'xe^xF^dx^ . (3) 
It is clear that under the parity transformation in three-dimensions 

x — > x 

y -»• -y 
t^t , 

the 9 term changes sign i.e 9 — ► —9. Therefore, the 9 term violates parity unless the free 
energy is an even function of 9. The field is a neutral field and hence the above term is 
gauge invariant. A simple integration by parts gives the following terms 

dxie^xF^x)^)) - (^f J- uxd\F^ u ( K x))(f)(x) . (4) 
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The dual of the field strength is defined as 

F x( x ) = 7f\^F^ v (5) 

The second term is seen to be 2d x F x (x)<ft(x) and this term can be non-zero in the presence 
of magnetic monopoles. In the presence of magnetic monopoles 

d x F x = m(x) (6) 

where rh(x) is the magnetic monopole density at x.The second term leads to the coupling 

26im{x)<f)(x) . (7) 

Another non-trivial contribution to Eq. ^| can come from a vortex line. This is seen by 
doing the integration by parts in Eq. ^| differently as 

2e fluX d ll (A u d X (j)) - 2e lluX A u d ll d x <f) (8) 

The second term can be written as 

- 29% J m^x)A^x) (9) 

The quantity 

m^x) = e^ Xu d x d u <p (10) 

is non-zero around a vortex line. Hence, the term in Eq. ^| is non-zero in the presence of 
magnetic monopoles and vortex lines and introduces new couplings between the topological 
( Thfj,(~kx) and m(*x) ) and the spin wave (A^ and <fi) degrees of freedom. In the absence 
of these topological excitations, the term in Eq. |^ will have no physical effect. Another way 
of motivating the term in Eq. |3| is by dimensionally reducing the four-dimensional term 
9 term in Eq. |2|. At high temperatures, the leading order contribution from a term like 
Eq. ^| is given by Eq. [3]. In this paper, we will present an analysis of the 9 term in three 
Euclidean dimensions. In three dimensions, the physical properties of this model are quite 
different from the four dimensional one. We will show that many of the interesting features 
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pointed out in |7], |] are also present in three dimensions. However, there are also significant 
differences. As already explained before, the 9 term becomes important whenever there are 
magnetic monopoles or vortex lines. On the lattice we can naturally define models which 
contain monopoles and vortex lines as possible excitations. This is possible if the degrees 
of freedom are considered as angular variables. We first present two well known lattice 
models which have monopoles and vortex excitations and then study the effect of the 9 
parameter in these models. We will see that the 9 parameter couples these two models 
in a non-trivial way and leads to a rich phase structure. The analysis presented here is 
based on the technique of duality transformations as applied to statistical systems. These 
techniques have been used very effectively to understand systems like the planar model ffl 



and the U(l) lattice gauge theory ||10|| . Recently in [[LlJ dual transformations have also 
been constructed for non-abelian lattice gauge theories. 

The main aim of this paper is to present an analysis of the 9 term in a lattice model 
containing both monopoles and vortex lines. The 9 term is introduced on the lattice by 
coupling the monopoles to the scalar field and the vortex lines to the gauge field. Many of 
the techniques used for studying the lattice model are well known in the literature but we 
present many details for the sake of completeness. Also, at some places we have managed to 
give more illuminating derivations of some of the steps in the analysis. The organization of 
this paper is as follows. In Sec. 1 and Sec. 2 we discuss lattice models which have monopoles 
and vortices as possible excitations. In Sec. 3 we present a detailed analysis of the model 
obtained by adding a 9 term. In Sec. 4 we make some concluding remarks. Some technical 
details are presented in the appendix. 

2 A Lattice model with monopoles 

The model considered here is given by the following action 

Si = —^^{d^(f) v - d v (f)^-2ir s^ u (x)) 2 + ip^m fl (x)4> fi (x) . (11) 

XpuV Xfl 



The partition function is given by 

Z 1 = trexp(S 1 ) . (12) 

The symbol tr denotes 



E E / ■ (is) 

m M =— oo s M „=— oo 00 

The variables <f)^ and m M are defined on the links of the lattice and the integer valued 
variables are defined on the plaquettes of the lattice. The symbol x denotes a three- 
dimensional vector and the symbol is the lattice derivative. The fields <fi^ are the 
gauge degrees of freedom whereas the integer valued variables s^ u are the monopole degrees 
of freedom. The above model describes gauge fields coupled minimally to current loops 
with the additional presence of magnetic monopoles. This model has the following gauge 
invariance 

^(aO-^ + ^A . (14) 

This also requires that 

<W = (15) 

Hence the summation over m M is restricted to a summation over closed loops. To see that 
this model describes magnetic monopoles consider the quantity 

F^(x) = d^ v {x) - d u (p^(x) - 2-ks^ v (x) . (16) 

Consider a configuration in which s fiu (x) 7^ on all the plaquettes pierced by the string 
in Fig. []. This figure shows the projection of the string on the XY plane (as broken lines) 
and the projection of the plaquettes ( which are in the YZ plane) on the XY plane (as solid 
lines). The string is on the dual lattice and begins at P and extends to infinity in the z 
direction. It is easily seen that for this configuration 

d\F\ = -2vr5(x P ) (17) 

which means that there is a magnetic monopole at P with a Dirac string in the z direction. 
Thus, the integer valued plaquette degrees of freedom account for magnetic monopoles. The 
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monopole and its associated Dirac string reside on the dual lattice. The monopoles in this 
model can be made explicit by means of a duality transformation. A duality transformation 
generally involves three steps. First, the quadratic part of the action is linearized by 
introducing an auxiliary field. The next step is to integrate over the original ( in this case 
</> M ) degrees of freedom and this will lead to a constraint for This constraint is solved 
by introducing degrees of freedom which are the dual variables and this solution is inserted 
back into the partion function. A few manipulations after that lead to the following form 
of the partition function 

Z\ = Tr exp ( — 8n 2 p g fh{x)G{x — x')rh(x') J exp j — ^- m^{x)G{x — x')ra^{x 

\ x,x> J \ x,x> 

X exp 2iri ^2 ra(x)G(x — x')e lll/ \d\K i 



The monopole density m(*x) is defined as 

fh{*x) = ^e\fivdxs^ u (x) . (19) 

The symbol *x denotes that the monopole is defined on the dual lattice point. This 
definition requires that the total monopole number in the system is zero. The symbol Tr 
denotes 

E (20) 

m M ,m=— oo,oo 

with the understanding that the summation is over configurations with zero total monopole 
number and closed loops of currents. The quantity K* u is a particular solution of 

d,Kl u =pm u . (21) 

In the above expression for the partition function G(x — x') is the three-dimensional Green's 
function which satisfies 

- d 2 G{x - x') = 5{x - x') . (22) 

The symbol d 2 is the lattice laplacian. The partition function in Eq. |l^ describes a gas of 
current loops (m M (x)) and magnetic monopoles (m(*x)). The monopoles and current loops 



interact among themselves with a three-dimensional Coulomb interaction. The last term 



in Eq. [18] describes the interaction a of a monopole with a current loop. This interaction 
is just the solid angle subtended by the current loop at the monopole and this is shown 
in the appendix. Before proceeding we briefly outline the steps leading to Eq. [H| The 
linearization of the action is accomplished by introducing an auxiliary field K^ v and the 
partition function becomes 



/OO J 
dK^trexjp(i ^ K^(x)(d^(f) u ~d u (f)^-2iT s^ u (x)) exp(— K^(x)) exp(ipJ2 "vO)</vO)) 

"°° X ilV ^Pg X XfMU 

(23) 



'2/3, 

Integration over 0^ results in the constraint 

2d tt K lw = pm v . (24) 

The solution of this constraint is 

2K^ = e^ x d x <j) + K; u . (25) 

Since 2K ilv has to be an integer the fields <f)(x) and K* u (x) are integer valued. However, 
the solution to the constraint is not unique because 

<f>(x) -> <f>(x) + c (26) 

also solves the constraint. Using this "dual" gauge invariance the range of integration over 
<p can be taken from — oo to oo. Substituting the solution of the constraint in Eq. ^5] in 
Eq. |23| and doing a gaussian integration over leads to Eq. |TJ. 

3 A Lattice model with vortices 

A lattice model describing vortices is defined by the action 

s 2 = 4 E(W-HWf + + i pE m W) • (27) 

" x [i r 
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The partition function is given by 

Z 2 = trexp(S 2 ) . (28) 

The symbol tr denotes 

£ / M ■ (29) 

m,( M =— 00 

In the above model the variables 9 and m are defined on the sites of the lattice and 
the variables n M are defined on the links of the lattice. 9 are the spin-wave degrees of 
freedom and are the vortex degrees of freedom. The above model describes vortex lines 
interacting with charged m(x) variables. The model has the following global invariance 

9{x) -> 9{x) + c (30) 

where c is any constant. This automatically requires 

J2m(x) = . (31) 

X 

The vortex lines in this model can be identified by considering the quantity 

Vp{x) = d^9{x) - 2nl^x) . (32) 

Consider a configuration in which / M is non-zero on all the links pierced by the world line of 
the string. This configuration is also shown in Fig. |6|. This configuration is a vortex running 
in the z direction which is the direction out of the plane of the figure. The accompanying 
string is chosen in the X direction and is indicated by the broken line which pierces the 
bonds between pairs of nearest neighbour sites (shown as solid lines). Any closed loop 
about this vortex line will give a non-zero value for 

£^(*) • (33) 

X 

The vortices in this model can again be made explicit by means of a dual transforma- 
tion. Introducing an auxiliary field as was done in the previous section and repeating the 
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procedure described before, we get an expression for the partition function 

p 2 

Wh 



Z 2 = Tr exp I —7^5- E m ( x )G(x — x')m(x') ] exp I —27i 2 [3 h m ll (x)G(x — x')m^{x r ) ] (34) 



x,x' 



X exp -2-ki m x (x)G(x - x^ex^d^Kl 



x,x' 



The Tr denotes 



£ (35) 

m,m=— 00 

The vorticity m^{*x) is defined as 

rh^-kx) = t^xdu^xix) . (36) 

From the above equation it is clear that the vortices form closed loops because 

dfjh^-kx) = . (37) 

As in the case of the monopole, the vortex lines and their associated strings live on the 
dual lattice. A closed vortex loop in the \xv plane will have a sheet swept by its string and 
the plaquettes in this sheet will be dual to the links with / M 7^ 0. The steps leading to the 
dual transformation are analogous to those in the previous section; only the auxiliary field 
is now introduced on every link. The constraint equation that has to be solved is 

d^K^x) = pm(x) . (38) 

The solution of this constraint equation is 

K^x) = e^dufa + K;{x) , (39) 

with K*(x) being a solution of the inhomogeneous equation. The dual gauge invariance in 
this case is 

4>(x) -> 4>(x) + d\A . (40) 

The first two terms in Eq. |34] describe the Coulomb interaction between m(x) variables and 
the vortex currents fh^-kx). {G(x — x') is the same Green's function as in Eq. ^2|) The last 
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term represents an interaction between the vortex currents and the m charges. K*(x) is 
the solution of the inhomogeneous equation 

d,K;(x) = pm(x) (41) 

This interaction is again proportional to the solid angle (apart from a negative sign) of the 
vortex current subtended at m. A demonstration of this fact is there in the appendix. 

Before we proceed to the model with a 9 term we can already see that the two models 
described above have a very similar structure. For instance, the monopole model has 

= (42) 

as the current conservation equation and 

J2 m(*x) = (43) 

X* 

as the monopole conservation equation. In the vortex model, the "m" charge conservation 
equation is 

]Tm(x) = (44) 

X 

whereas the vortex conservation equation is 

d^m^-kx) . (45) 

The roles of the conservation of ordinary charges and topological charges are clearly reversed 
and the two models are dual to each other. This duality will be made more precise in the 
next section. 
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4 The Coupled Model 

Now we can couple the two previous models by introducing a 6 term as explained in the 
introduction. The 6 term is defined by introducing two additional couplings as 

S e = 7pE™M#*0 + {-12™»(* x )M x ) ( 46 ) 

' X,fM X,fM 

The full action of the coupled model is given by 

S = S 1 + S 2 + S e , (47) 

which can be written out as 

S = Z lrY, ( 5 A ~ " ^s, u {x)f - ^ J2(d,9(x) - ^n,{x)f + (48) 

x \xv x n 

ip^m^x)(j) v {r : ) + ipJ2m(x)9(x) + |^ rh{-kx)4){x) + rh^-kx^^x) 

X jl X ^" X,[l ^" X,[l 

The partition function of this model is given by 

Z = trexp(-S) . (49) 
The trace represents the following sum over states: 

00 00 />00 PCO 

E E / / • (so) 

/ \ / \ J —CO J —CO 

m(x) —com^x) -co 

When 9 = the two models are decoupled and the partition function is simply a product 
of their separate partition functions. 

Z = Z X Z 2 . (51) 

The model at 6 = represents a system of monopoles, currents, m-charges and vortices. 
However, the excitations in one system do not interact with those in the other system. 
The model at 6 = is seen to be trivially self-dual. This follows by noting that the 
transformations 

* i^k (52) 
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yn 16n 2 p g 
m(x) — > rh(*x) 

m^x) -> rh^-kx) 



leave the partition function unchanged. The dual transformation maps every point on the 
hyperbola 

AA = ^ (53) 
onto itself. The region (3 g (3h < is mapped onto the region f3 g Ph > by the dual 
transformation and vice-versa. However, this self duality property is trivial because the 
system on which it acts is a product of two decoupled systems. Nonetheless we mention 
this fact here because, as we will show later, at certain values of 9 the self-duality will still 
persist When 9^0 the two systems are coupled in a non-trivial way. There is a cross 
coupling between the spin-wave excitaions of one syatem and the topological excitations 
of the other system. For instance, the monopole rh(*x) couples to the spin-wave field 
m(x) and similarly the vortex m^x) couples to the gauge field 0^. The main point of this 
paper is that this coupled system defines an interesting model that has some exact duality 
symmetries. It also has a rich phase structure as a function of 9. The analysis previously 
described for the monopole and vortex models can be repeated in the same way as before. 
The only point to note is that the monopoles and vortices are defined on the dual lattice 
whereas the gauge fields and the spin variables are defined on the original lattice. However, 
we can approximately take the point on the dual lattice to coincide with the point on the 
original lattice. There will be corrections to this approximation but these will involve higher 
derivative terms which can only affect the short wavelength behaviour of the system. The 
9 term changes the constraint equation for the two auxiliary fields. The new constraint 
equations become 

p9 _ 

d^K^ = pm u + —m^-kx) (54) 
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duKJx) = pm(x) + —fh{*x) . 

lix 

The only change is in the inhomogeneous part of these equations 

p9 „ 

d n K U = P m » + ^ m ^ x ) (55) 
p9 

duK*(x) = pm(x) + —m(*x) . 

2,7V 

Repeating the steps performed for the monopole or the vortex model leads to a 6 dependent 
partition function 



Z e = tr(56) 



(p 2 9 6 s 

~gp- g Y,( m »( x ) + 2^™»( X )) G ( X - x ') i m »( x ') + 2^™»( x> ) 

exp f — y2(m(x) + —m(x))G(x — x')(m(x') + — ra(x') ) exp f — 8ir 2 f3 g y^ j m(x)G(x — x')m{x 
exp I -27r 2 /3 h ^ m il {x)G{x - x'jrh^x 1 j exp 2m ^ m(x)G(x - x')e iJlv xd\K i 

\ xx' / \ x,x' 

exp — 27TZ m x (x)G(x - x^ex^d^K* 

\ x,x' 

It can be seen that the 9 term couples the gauge and spin models and introduces additional 
interactions in each of them. The partition function of the model is no longer separable 
into a spin part and a gauge part and cannot be written as 

Z e ^ Z,Z 2 (57) 

Hence, the phase structure of this model can be quite complicated. One of the immediate 
consequences of the term is that the term representing the Coulomb interaction between 
electric loops gets modified. This means that the vortex loops (for which fn^-kx) ^ 0) 
acquire an electric charge which is given by 

Q 

Qn(x) = m^x) + — fn^-kx) . (58) 

For instance, a loop having only vorticity (m^x) = 0, m^x) ^ 0) will still have an 
electric charge given by Jj^fh^-kx). Likewise the Coulomb interaction between m charges 
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is also modified and acquires a piece due to the monopoles 

Q 

Q(x) = mix) + — m(*x) . (59) 

Z7T 

Similarly a point having only non-zero monopole density will have an additional "m" charge 
|^m(*x). The electric charges of the vortex loops and the "m" charges of the monopoles 
can now take fractional values because of 9. It is convenient to associate a vorticity and 
an electric charge to every closed loop. The values of these charges are plotted in Fig. ^| 
The first thing to notice about this model is that the partition function is periodic in 9. 
This follows from the simple fact that we can always shift the summed variables and 
m as their summation range is infinite. In the presence of a 9 term the model is no longer 
self-dual under the transformations in Eq. [52[ However, for certain specific values of 9 the 
model is still self dual. These are the values for which 2ir/9 is some integer q. To see the 
self duality for these values of 9, we make a simple change of variables 

mv + = JjL (60 ) 
q q 

m I 

m H = — . 

q q 

Expressing the partition function in terms of l u ,m u , I and m, the partition function reduces 
to the the original one provided we identify 1^ with m M , I with m, and make the following 
changes. 

-» ^ 
^ 9 i6ttV/^ 

p 2 

^ ~* I67r 2 q 2 g ' 

Unlike the dual transformations in Eq. [52| the above dual transformations contain more 
information because they now act on a system which is not a decoupled system. These 
dual transformations can be used to understand the phase diagram of the model in one 
region if the phase diagram is known in another region. The points of the phase diagram 
which are left invariant under the dual transformation are the points on the hyperbola 

2 

\§-K l q l 
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The self-duality in Eq. |6l] holds true only when 9 = 2rc/q,q being an integer. We now 
point out another symmetry that is present in the model for arbitrary values of 9. This 
symmetry can be deduced by requiring the partition function to be invariant under the 
following transformations: 

m — > — m (62) 
9^9' 

As the partition function is a trace over the m, m, m^, degrees of freedom it can be 
written as 

Z = Zg(m, m, m M , m M ) (63) 

Imposing the condition 

Z e t(-rh, -m, rhp, m v ) = Z e (m, m, rh^) (64) 
we get the following set of equations 

"i+Sk-i (65) 

P2 W + p26 ' 



2(3 h ^ Sir* ft 
2 p 2 9 2 f_ 
Ph + 32n 2 f3 g ~ 8/3' g 

8/3 9 " 32^7 + Pk 

We get four equations in the three variables, P g ,Ph, and 9. Interestingly, these four equa- 
tions are consistent and have the following solution: 



_ 4:TT 2 (3 h p 2 
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Pg Q4Tr*{3 g f3 h +p 2 9 2 
2 _ 16ttW 

(647r 4 /3 9 ^+P 2 ^) 2 

This means that we can define an exact duality symmetry in the model which now acts 
on a three coupling space. The duality equations can be recast in a more concise form by 
defining the variable 

z = yfWh) ■ (67) 
In terms of the variables z and 9, the duality equations are 

Z (64tt 4 z 2 + p 2 9 2 ) 2 { ' 

2 16n 4 p 4 2 



0' 

(64tt 4 2 2 + p 2 6 2 ) 2 
If we define the complex coupling ( by 



p 2n 

the duality transformation can be expressed as 



4irz id , . 

C = + 77- , (69) 



C^i • (TO) 
Since the action is also periodic in 9 with period 2ir, the transformation 

C^C + * (71) 

is also a symmetry of the model. These two transformations do not commute with each 
other and generate the group SL(2, Z). This symmetry group was first pointed out in [[§] in 
the four dimensional model. In the three-dimensional model that we are considering here, 
the duality symmetry has a slightly more complicated form as in Eq. ^6] and the group 
SL(2, Z) acts on the variable z which is given by Eq. [67]. 

We now proceed to study the phase structure of this model. The phase structure of 
this model will be a function of the three parameters (3 g ,(3h and 9 at a given value of 
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p. The various phases of the model will be characterized by the behaviour of the the 
m(x),rh(x),m ll (x),rh fl (x) excitations. Depending on the density of these excitations, the 
phases in this model will have different physical properties. In order to arrive at the phase 
structure of the model we will use simple free energy arguments based on the energy and 
the entropy of loops. Though these arguments are admittedly crude and ought to be 
substantiated by other methods they provide us with a qualitative picture of the phase 
structure. We first note that there will be regions in parameter space where the loops will 
condense. There are two kinds of loops that can condense. They are labelled by m M (r) 
and m^(r); m M (r) loops are the world lines of electrically charged particles, m M (r) loops are 
the vortex excitations. A composite loop having both m M (r) and m M (r) is also possible 
and it will be referred to as a composite loop. These composite loops are formed by the 
binding of a vortex to an elecrtric current loop. In addition to these loops there are also 
the point like excitations labelled by m(x) and rh(x). rh(x)s are the magnetic monopoles 
and m(x)s will be referred to as spin charges. The reason for this terminology is that 
the m(x)s represent source terms for the spin variables. In addition to these excitations 
there are also excitations which simultaneously carry magnetic charge and a spin charge. 
These composite objects will be referred to as composite charges. The composite charges 
are formed by the binding of a magnetic monopole to a spin charge. A composite loop or 
a composite charge is a combination of an ordinary charge and a dual excitation. If we 
neglect the long range Coulomb interaction, a crude estimate for the free energy of a loop 
of length L having charges (m^,rh^) is 

2ttV 

F(L) = (^(m+^rh) 2 + 2n 2 f3 h (m) 2 )G(0)L-(\ogc)L . (72) 

In the above approximation, only the self interaction of the loop contributes to its energy. 
Condensation of loops having charges (m M ,m M ) occurs whenever the loop free energy be- 
comes negative. The constant c depends on the co-ordination number of the lattice and is 
approximately 5 for a three-dimensional cubic lattice. In the above approximation only the 
self energy of the loop has been considered in the expression for the free energy. Though 
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this is quite a severe approximation, we expect it to reproduce the general features of the 
phase diagram. As we have already noted, the effect of the 9 term is to give an electric 
charge to the vortex lines as 

Q 

q u = m v + —rh v . (73) 
Ztt 

Similarly, the magnetic monopoles get an effective spin charge 

e . s 

q — m-\ m . (74) 

Itx 

We first consider the phase diagram at 9 — 0. When 9 = the partition function can be 
decoupled into two partition functions 

Z = Z 1 Z 2 (75) 

and each of these systems can be considered individually. If we consider Z\ first, the free 
energy of a loop of length L carrying electric charge m M is given by 

2ttV 

F(L) = —f-m 2 G{$)L - (log c)L . (76) 

Pg 



Thus F(L) < if 



P > 2nVG{0) (77) 
' 9 logc ' 1 ] 



Therefore the model described by Z\ will exist in two phases, a small f3 g phase in which 
the current loops are very sparse and a large f3 g phase in which the current loops are very 
dense. The other excitations present in this model, namely the magnetic monopoles will 
have a density given by 

p((3 g ) = exp(-87r%G(0)m 2 (x)) . (78) 

Unlike the density of current loops, the magnetic monopole density falls continuously to 
zero as the coupling constant f3 g is increased. We can make a similar analysis of the model 
Z 2 . The same free energy arguments applied to the vortex loops gives F(L) < if 

log(c) 

^27^(0) • 1 j 
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The small fih phase has a high density of vortex loops whereas the large /3h phase has a 

very low density of vortex loops. The spin charges m change continuously with fih as 

, p 2 G(0)m\ . . 

p(p h ) = exp(- p V ) • (80) 

For the coupled model the analysis can be repeated but there are now three coupling 



constants (3 g ,Ph and 9. We consider the condensation condition given by Eq. [72] for various 
limiting values of (3 g and fih- It we fix 6 = tt the different excitations present in the system 
are those which correspond to the black points in the charge lattice in Fig. |6|. The excitations 
corresponding to these black points are multiples of these fundamental excitations. 

1. Electric current loops (1,0): These have a non-zero electric charge. 

2. Vortex loops (0, 1): These have a non-zero electric charge because of Eq. [731 . 

3. Composite loops (1, —2): These loops have exactly zero electric charge again because of 



A similar charge lattice can be drawn for the point like excitations in the model. The 
excitations of this charge lattice are multiples of the following fundamental charges. 
l.Spin charges (1,0): These have only "m" charge. 

2. Magnetic monopoles (0, 1): These have a non-zero "m" charge because of Eq. |74|. 

3. Compositely charged point like objects (1, —2): These have a zero "m" charge because of 
Eq. |7|. 

We have to compare the free energies of these condensates and then choose the one with 
the lowest free energy. First we can explore various limits of the coupling space. 
1-Ph = o. 

The free energies of the possible condensates are given by 

F(L)/L = ^G(0)(m+im) 2 -log(c) . (81) 

(a) (l, ):^G(0)-log(c) 

(b) (0, 1) : ^G(O) - log(c) 
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(c) (1,-2) :-log(c) 

The above estimates show that the composite loops (1, —2) will always have the lowest free 
energy. Therefore, in this limit, there is always a condensate of composite loops (1, —2) 
and there will be no phase transition on this axis. 

2. j3h = oo The general free energy relation Eq. [T2| shows that this limit forces « 0. The 
free energy condition becomes 

F(L) = |U 2 G(0) - log(c) . (82) 

Since fh is forced to be zero we only have to consider electric (1, 0) loops. F(L) of (1, 0) 
loops becomes negative for f3 g > g^j^y- This implies that there is a phase transition on this 
axis from a small (3 g phase containing very few current loops to a large (3 g phase containing 
very large current loops. 

3. /3 g = 0. In this limit we get the following constraint 

(m + ]-m) = . (83) 

This leaves only the composite loops (1, —2) as possible excitations. The free energy con- 
dition becomes 

F(L) = 8n 2 p h G(0) - log(c) . (84) 
Thus, (1, —2) loops will condense for 

log(c) 

Ph< 8^G{0) ■ (85) 
Therefore, on this axis we expect a small (3h phase in which (1, —2) loops condense and a 
large j3h phase in which the density of (1,-2) loops is very small. 

4. /3 g = oo. The free energy condition is 

F(L) = 27T 2 (3 h m 2 G(0) - log(c) . (86) 

(a) (l,0):F(L) = -log(c) 

(b) (0, 1) : F(L) = 2vr 2 /3,G(0) - log(c) 
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(c) (l,-2):F(L) = 87r 2 /3,G(0)-log(c) 

It is clear that the condensate (1, 0) always has the lowest free energy. Therefore, the 
electric current loops are always dense on this axis and there is no phase transition on this 
axis. The other special case to be considered is 

5 A = p h 

In this limit the free energy condition becomes 

F(L) = f^(m + im) 2 G(0) + 2n 2 (3m 2 G(0) - log(c) . (87) 
op 2 

The condensation criterion can be written as the interior of the ellipse on the charge lattice 

2 ~ 2 

,q m\ , , 

& + -jT><l • (88) 

The major and minor axes of the ellipse are given by 

8/?log( C ) 

G G(0)p 2 m 

2 = log(c) 

2vr 2 /?G(0) 



When a » 1 and b << 1 the ellipse is very flat along the q-axis and the condensate 
with the lowest free energy is (1,0). When a << 1 and b » 1 the ellipse is very flat 
along the rh axis and the condensate with the lowest free energy is (1,-2). For values 
of a and b such that they are comparable, the free energy is the lowest for the (0, 1) 
condensate. Therefore we generically expect three phase transitions on this axis. This case 
occurs in J7| in the four- dimensional model. It is interesting to note that exactly the same 
condensation condition appears in the three-dimensional model with the only difference 
being that the four- dimensional Green's function gets replaced by the three-dimensional 
one and the entropy factor takes different value. 

Now that we have discussed the condensates in the various regions of parameter space, 
we can propose the phase diagram of this model at 9 = tt. It is given in Fig. || So far we 
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have only discussed phase transitions of the electric, vortex or composite loops. We have 
already noted that there are other excitations present in this model. These are the point 
like excitations namely the magnetic monopoles, "m" charges and the composite charges 
of m and m.The density of these excitations changes continuously with (3 g ,(3h and 9. The 
density of the magnetic monopoles and the "m" charges are given by 

-v 2 1 

p m> rn = exp(-87r 2 ^G(0)m 2 ) exp(-f-G(0)(m + -mf) . (90) 

This is one of the special features of the model in three-dimensions. As these excitations 
are point like, they form a three-dimensional Coulomb gas of charged objects whose density 
is always non-zero. The densities of the three possible types of point like excitations in this 
model are give by 



1.(0, 1): Magnetic monopoles 



2.(1,0): Spin charges 



2 

p^exp-(87r 2 /? g + |-)G(0) . (91) 



2 

Pm « exp(-^-G(0)) (92) 
^Ph 



3.(1,-2): Composite charges 



p c w exp(-647r 2 /3 ff G(0)) (93) 



Since each of these densities is always non-zero, we have a three- dimensional Coulomb gas 
containing three species of charged particles. The densities of these three species of charged 
particles determine various correlation functions (which are defined later on) independently 
and can lead to different effects at different values in the parameter space. The behaviour 
of physical correlation functions will be determined by the density of electric and vortex 
loops as well as the density of magnetic monopoles and "m" charges. The various regions 
in the schematic phase diagram in Fig. || are discussed below. 

A. In this region only the composite loops (1, —2) condense. At large values of (3h, the 
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density of magnetic monopoles and spin charges is large small whereas at small values of 
f3h they become non-zero but obey the relation 



The density of magnetic monopoles decreases as (3 g is increased. 

B. In this region the electric current loops (1, 0), the vortex loops (1, 0), and the composite 
loops (1, —2) have almost zero density. Therefore, this phase is essentially free of all loop 
excitations. The magnetic monopoles (0, 1) and spin charges (1,0) have a large density at 
the top left corner of the phase diagram which decreases as as we go away from it in any 
direction. 

C. In this region the electric loops (1, 0) condense whereas vortex loops (0, 1) have very low 
density. The magnetic monopoles (0, 1) have a very low density and the spin charges (1, 0) 
have a large density at the top tight hand corner of the phase diagram. 

D. In this region the vortex loops (0, 1) condense and the electric current loops (1, 0) have 
a very low density. The magnetic monopoles (0, 1) and spin charges (1, 0) both have a 
density which is not very large or very small. 

We will now show that these different phases can be characterized by the behaviour 
of correlation functions which are of the order disorder type. Before we do this we briefly 
describe how these correlation functions are defined in each of the models discussed in the 
previous two sections. The correlation functions are simple generalizations of the Wilson 
loop |I2"|| . In the model described in Eq. [18] we can introduce an external current J^x) 
in a loop C on the lattice and a monopole-antimonopole pair at points *X\ and *X2- The 
monopole pair is introduced by choosing = 1 on a string joining *x\ and -kX2- The 
correlation function is defined by 



m H = 

2 



(94) 



W(C,x h x 2 ) 




V + s'(x)) 



(95) 



This is equivalent to making the following change in m(*x) 



rh(*x) 



rh(*x) + p(*x) 



(96) 
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where p is defined as 

p{*x) — 8{kX — — 5(*X — *X 2 ) ■ (97) 

The correlation function W measures the free energy of an external current loop J M on 
C and an external monopole-antimonopole pair at points -kxi and *x 2 - We can similarly 
define a correlation function for the model described in Eq. [34] as 

Z 2 (m^m + M(x)J u ^L(x)+l'(x)) , . 

W(*C,x h x 2 ) = — ^ . (98) 

Here I' ^ on all the links which are dual to the plaquettes in a surface bounded by the 
loop -kC on the dual lattice and M(x) has the form 

M(x) = 6(x - x x ) - 6(x - x 2 ) . (99) 

The correlation function W measures the free energy of an external vortex loop *C and 
an external pair of spin charges at x\ and x 2 . The correlation functions W and W are 
examples of order disorder variables as they involve both the gauge(spin) degrees of freedom 
and the monopole (vortex) degrees of freedom. In the interacting model we can consider 
the following correlation function 



W int (C, *C, xx,x 2 , *x 2 ) 



Z(m M -» + Jft, -» + s'^ u , m -» m + M, -» l M + ^) 

(100) 



In the presence of interactions (when # 7^ 0) 

Wi n ^ WW (101) 

where we have suppressed the arguments of the correlation function. This correlation 
function obeys the same duality invariance as the partition function 

Wint{C,*C,Xi,X2,*Xi,*X 2 )l3 a ,l3 h ,0 = W int (*C, C, *X X , *X 2 , %1, X 2 ) p' g ,/3' h ,e' (102) 

where j3',f3' h and 0' are related to the unprimed values by the duality equations in Eq. 
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It is well known that a dilute gas of magnetic monopoles results in an area law behaviour 
for the Wilson loop . This effect arises from the form of the monopole-current loop inter- 
action. Also, an external monopole-antimonopole pair will experience a screened Coulomb 
potential in this gas. Since the form of the spin charge-vortex loop interaction is the same 
as the monopole-current loop interaction, a dilute gas of spin charges will result in a similar 
area law behaviour for an external vortex loop. Again, an external pair of m charges will 
experience a screened Coulomb interaction in this gas. Similarly, composite charges (1, —2) 
will result in an area law for an external composite Wilson loop that consists of a loop of 
electric charge (1) and vorticity (-2). However our model also has loop like excitations 
which can screen external Wilson loops and vorticity loops which can screen external vor- 
tex loops. For instance, a condensate of (1, 0) loops can screen a Wilson loop of charge 
(1,0) and yield a perimeter law. When both current loops and monopoles are present, the 
Wilson loop will have an area law piece but the perimeter law piece will dominate at large 
distances. Similarly, when (0, 1) vortex loops condense an external vortex loop of charge 
(0, 1) can be screened by these loops and again result in a perimeter law for the external 
vortex loop inspite of the the presence of m charges which alone would have yielded an 
area law. The same holds true for composite (1, —2) loops and (1, —2) charges. From the 
partition function in Eq. |5B] we see that there is no monopole- vortex or spin charge- current 
loop interaction. Thus these excitations cannot influence each other in any drastic way. 

5 Conclusions 

In this paper we showed that it is possible to define an analogue of the 9 term in three 
dimensions. We presented our analysis of the 9 term in a lattice model which was motivated 
by the work of the authors in f7|, || . The 9 term couples magnetic monopoles to the scalar 
field and vortices to the gauge field. The phase structure of the model changes as a function 
of 9. In fact, the interactions in this model arise entirely because of the non-zero value of 
9. This model has an exact duality symmetry which acts on a three parameter space. This 
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seems to be the first example of a statistical model in which the duality transformation acts 
on a three coupling space. We also made a qualitative analysis of the phase diagram using 
energy entropy arguments and showed that at non-zero values of 9 there are phases in which 
excitations having composite loops condense. A special feature of the three-dimensional 
model is that there are novel point like excitations of different species which form a three 
dimensional Coulomb gas. Our analysis is also instructive in understanding how the dual 
transformation works in systems containing both vortices and monopoles. Since the phase 
diagram of this model (which was studied at 9 = ir) admits phases with charged vortices, 
it will be interesting if some condensed matter systems can be described by this model. 
Another interesting avenue is to introduce 9 like terms in non-abelian gauge theories and 
look for oblique phases at non-zero values of 9. 
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6 Appendix A 

In this appendix we will examine the form of the monopole- current loop and the vortex-m 
charge interaction. The monopole-current loop interaction is given by 

J2 rh(*x)G(x - x')e^K uX {x') (103) 

x,x' 

K* is a particular solution of 

d,K^ = pm u (104) 

Consider a configuration in which the monopole is along the z — axis and at a distance z 
and there is a circular current loop (of radius R) is in the X-Y plane i.e only m x ,m y ^ 0. 
It is easily seen that a particular solution of Eq. |104| is 



K xy = —p insideC (105) 



K xy = otherwise 



K xz = Kyz = 



This reduces the monopole-current loop interaction to 

- 2 rh{*x)d z ,G{x - x')K* y (x') (106) 

x,x' 

The co-ordinates of x and x' are given by 

x (0,0, z) (107) 
x' (x',y',0) 



Using 



OrCU -.,■') = \ X _ ^ (108) 



\x — x 



the interaction is given by the expression 



r z dx'dy' 
2pm J. (x» + y» + z»)W (109) 
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The integration is over the area of the loop. The evaluation of this integral is straightforward 
and gives 

/ = 2pm(2n(l - cos(a)) (110) 

where a is the azimuthal subtended by the current loop at the monopole. 
Now we turn to the vortex-m interaction. It has the form 

-J2rh x (x)G(x - x')e x ^Kt{x') (111) 

x,x' 

K* is the solution of the inhomogeneous equation 

d„K*{x) =pm(x) (112) 

Consider a configuration in which m (x) 7^ at a point on the z axis a distance L from the 
origin. The solution of the inhomogeneous equation is chosen to be 

K = k; = o (us) 

K* z = P e(z-L) . 



The G function has the property 



Q{z) = 1 z > (114) 
e(z) = z<0 . 



The co-ordinates of the points x and x' are given by 



x (Rcos(9),Rsm(9),0) (115) 
x' (0,0, z') 



Again using Eq. |108| we are led to the integral 

- 2 ™P R2 J L {R 2 + 3,2)3/2 ( 116 ) 

This integral is also straightforward to evaluate and gives the solid angle interaction 

I = -2pm(2n(l - Cos(a)) (117) 
29 



which is just the negative of the monopole- current loop interaction. 

The purpose of this appendix was to show that the monopole- current loop interaction has 
the same form as the interaction between the m charge and the vortex loop apart from a 
sign factor. 
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Figure 1: The string attached to the dual excitations runs on the dual lattice. 



Figure 2: Charge lattice at 6 = n. 



Figure 3: Schematic phase diagram at 6 = 7r.The X-axis is the (3 g coupling and the Y- 
is the Ph coupling. 



32 



yt 



A 



c 



fa t» 



a 



q 



